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■  For  diffuse  targets,  features  such  as  shape,  size,  and  motion  can  be 
determined  from  a  time  series  of  images  from  either  angle-angle  passive 
telescopes  or  range-Doppler  radars.  The  extracted  target  features  can  then  be 
used  for  automated  target  recognition  and  identification. 

An  algorithm  that  uses  scene-analysis  techniques  has  been  developed  to 
perform  the  feature  extraction.  The  algorithm  first  processes  the  images  to 
suppress  noise,  then  applies  a  two-dimensional  slope  operation  for  edge 
detection  to  determine  the  target  boundaries.  Next,  Hough  transforms  are  used 
on  the  target  edges  to  detect  straight  lines  and  curves,  which  are  subsequently 
refined  with  line  and  curve  fits.  Groups  of  the  fitted  lines  are  then  examined  to 
form  cylinders  and  cones  representing  typical  target  components.  After  these 
shapes  have  been  identified,  the  target  configuration,  size,  location,  and  attitude 
can  be  estimated.  The  target  motion  can  then  be  inferred  from  a  time  series  of 
attitudes  that  have  been  extracted  from  a  sequence  of  images. 


For  a  target  with  rough  surfaces,  electromag¬ 
netic  signals  are  reflected  and  returned  from 
scatterers  that  are  distributed  over  the  entire 
target  surface.  The  resulting  imagery,  whether  angle- 
angle  (passive  telescopes)  or  range-Doppler  (radar), 
will  show  a  diffuse  object  with  surface  returns  that 
became  apparent  along  the  sensor  line  of  sight  (LOS). 
From  such  images,  the  target  shape,  size,  and  orienta¬ 
tion  can  be  determined  from  pattern  recognition  and 
identification  techniques.  And,  with  a  time  sequence 
of  images,  the  motion  of  the  target  can  be  estimated 
from  its  orientation  history. 

In  the  feature-extraction  algorithm  developed  at 
Lincoln  Laboratory,  the  images  are  first  processed  to 
s.jppress  noise  and  to  smooth  the  image  surface.  Edge 
detection  is  then  performed  to  determine  the  target 
boundaries,  and  the  detected  edge  points  are  inte¬ 
grated  into  line  segments  to  form  target  shapes.  Next, 
target  dimensions  and  orientations  are  measured  from 


the  line  representation  of  the  target,  which  is  assumed 
to  be  axisymmetric. 

In  this  article,  examples  of  target  feature  extraction 
are  demonstrated  for  both  angle-angle  and  range- 
Doppler  images.  For  angle-angle  images,  target  orien¬ 
tation  is  obtained  from  the  projected  elliptical  shape 
of  circular  components  and  the  projected  length  of  a 
symmetric  body  axis.  For  range-Doppler  images,  tar¬ 
get  dimensions  are  first  determined  from  a  sequence 
of  range  and  Doppler  extents  that  have  been  extracted 
from  the  images.  The  target  aspect-angle  history  can 
then  be  derived  from  that  same  sequence  of  range  and 
Doppler  extents  (as  for  angle-angle  images)  by  using 
the  estimated  target  dimensions. 

Image  Scene  Analysis 

The  goal  of  target  feature  extraction  is  to  obtain  shape 
information.  From  such  information,  the  size,  orien¬ 
tation,  and  position  of  a  target  can  be  estimated. 
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FIGURE  1.  Feature-extraction  process. 


Shapes  are  composed  of  lines  and  curves  that  are 
collections  of  boundary  edge  points.  Such  edge  points 
can  be  located  with  edge-detection  techniques.  The 
detected  edge  points  can  then  be  associated  with  lines 
and  curves  by  using  predetermined  fitting  constraints. 
For  the  current  application,  cylindrical  and  conical 
shapes  are  considered,  and  Hough  transforms  are  used 
to  detect  the  presence  of  lines  and  curves  [1]. 

Figure  1  illustrates  the  feature-extraction  process. 
First,  image  processing  is  performed  to  enhance  the 
edge-detection  process.  The  image  processing  includes 
the  application  of  median  filters  to  remove  isolated 


noise,  the  averaging  of  multiple  image  frames  to  en¬ 
hance  the  image  signals,  and  the  spatial  smoothing  of 
the  image  surface.  Thresholds  are  then  applied  to 
remove  the  image  background.  Next,  a  linear  Hough 
transform  is  used  to  detect  and  collect  lines,  which  are 
later  examined  to  form  the  boundaries  of  cylinders 
and  cones.  The  position  and  orientation  of  the  axis 
of  symmetry  of  each  detected  shape  are  then  deter¬ 
mined.  For  angle-angle  imagery,  an  elliptical  Hough 
transform  is  generally  used  to  detect  and  fit  the  tar¬ 
get-base  curve  to  r.l’ow  for  the  subsequent  calculation 
of  the  target  aspect  angle.  For  range-Doppler  imagery. 
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an  elliptical  Hough  translorm  is  used  to  determine 
the  target  dimensions,  and,  hence,  to  obtain  the  tar¬ 
get  aspect  angle. 

1  he  three  tnain  steps  ol  scene  analysis  tor  the 
extraction  ot  target  features  are  edge  detection,  line 
identification,  and  curve  extraction.  These  three  steps 
are  described  in  the  following  subsections. 

Fdge  Detection 

1  he  boundaries  of  objects  in  an  image  exist  at  loca¬ 
tions  where  the  image  values  change  abruptly.  These 
abrupt  changes  can  be  detected  by  a  spatial  difference 
operator  [2].  I  he  design  ot  sueh  an  operator  depends 
closel)-  on  the  quality  and  complexity  of  the  image, 
and  on  the  desired  level  of  feature  extraction.  The 
current  application  uses  spatial  gradient  operators  [3]. 

The  operator  center  is  placed  at  an  image  location, 
and  the  result  'd  ilic  c  'nvolution  operation  on  the 
image  values  repre.sents  the  local  gradient  at  that  im¬ 
age  location  and  in  the  operator’s  direction.  The  edge 
value  at  the  image  location  /  (.v,  /)  can  be  expressed  as 

G[l{x,y)]=  +C;  , 

where  G\,  and  are  the  image  gradients  in  the  hori¬ 
zontal  and  vertical  directions,  respectively.  The  edge  is 
m  the  gradient  direction,  which  can  be  obtained  as 

B  =  tan”'  — , 

where  is  measured  with  respect  to  the  x-axis. 


Ciradient  operators  ttf  size  5,  as  shown  in  Table  1 , 
are  used.  In  general,  the  operator  size  should  be  de¬ 
creased  for  coarser  image  resolutions  because,  as  the 
resolution  becomes  coarser,  the  object  boundaries  be¬ 
come  closer  to  each  other. 

The  calculated  edge  values  are  small  for  a  smooth 
image  surface  and  large  for  a  discontinuous  surface. 
Thus  edge  points  can  be  detected  by  applying  a  thresh¬ 
old  to  the  edge  values.  The  threshold  can  be  varied, 
depending  on  the  level  of  detail  desired  for  edge 
extraction.  The  examples  presented  in  this  article  use 
a  threshold  equal  to  10%  of  the  dynamic  range  of  the 
edge  values,  i.e.,  G,,,;,,  +  0.10(6',,,,^  -  The 

choice  of  the  threshold,  however,  can  be  optimized 
with  respect  to  the  histogram  of  the  edge  values. 
Because  of  the  spatial  extent  of  the  operator,  the  edge 
values  obtained  after  the  application  of  a  threshold 
generally  are  thickly  populated  near  or  at  target 
discontinuities.  Thus  a  procedure  for  non-maxima 
suppression  can  be  applied  for  the  further  thinning  of 
these  edges  [4],  In  the  procedure,  the  edge  value  at  a 
point  is  set  to  zero  if  the  value  is  not  the  local  maxi¬ 
mum  in  the  direction  perpendicular  to  the  edge  direc¬ 
tion.  Figure  2  shows  an  example  of  the  detected  edge 
points  superimposed  over  a  gray-scaled  image.  The 
edge  points  include  the  outline  of  the  cones  and  the 
silhouette  of  the  cylinder  against  the  background. 

Line  Identification 

After  the  edge  points  have  been  detected,  they  must 
be  associated  with  lines  and  curves.  For  simple  closed 
contours,  the  edge  points  can  be  chain  coded  [5],  and 
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FIGURE  2.  Gray-scaled  image  with  detected  edge  points. 

the  shape  description  can  be  performed  with  syntactic 
pattern  grammars  [6].  The  current  application  con¬ 
siders  images  of  complex  objects  comprising  a  combi¬ 
nation  of  shapes,  from  which  broken  edges  and  edges 
inside  other  object  boundaries  are  permitted.  Thus  a 
layered  approach  for  shape  formation  is  required. 

Hough  transforms  can  be  used  to  detect  the  lines 
and  curves  in  an  image.  In  Figure  3(a),  a  line  in  two- 
dimensional  (2-D)  space  is  represented  by  the  di¬ 
rected  orthogonal  distance  h  to  the  origin  and  the 
angle  6  that  h  forms  with  the  x-axis.  Any  point  (x,  y) 
on  the  line  will  satisfy  the  following  equation: 

h  =  X  cos  &  +  y  sin  0 . 

The  set  of  lines  passing  through  a  given  point  (x,  y) 
can  be  plotted  in  h-i)  space.  The  result  will  be  a 
sinusoidal  curve,  as  shown  in  Figure  3(b')  Nore  rhar 


the  h  curves  ot  two  points  in  .v-y  space  w  ill  inierseci 
perit)dically  in  h-U  space,  and  the  points  ot  intersec- 
tit)n  (spaced  .t  radians  apart)  will  correspond  to  the 
line  defined  by  the  two  points  in  .v-y  space.  Fhus,  it 
the  /'curves  ot  all  edge  points  in  an  image  are  plotted 
in  />-<> space,  the  points  at  which  the  /'curses  intersect 
will  correspond  to  prominent  lines  in  the  image. 

In  practice,  the  />-^Hpace  is  divided  into  accumula¬ 
tor  cells,  and  the  cell  value  ot  a  particular  cell  is 
increased  by  1  each  time  a  curve  passes  the  cell's 
location.  Because  the  direction  ot  each  edge  point  is 
known  from  the  edge-detection  process  (discussed 
earlier),  onfy  the  portion  of  the  h  curve  tor  angles 
around  that  direction  needs  to  be  searched.  A  thresh¬ 
old  is  chosen  at  the  halt-length  (in  pixels)  ot  the 
shortest  line  expected  to  be  extracted  trom  the  image. 
Only  those  accumulator  cells  whose  cell  values  are 
greater  than  the  threshold  are  selected  tor  further 
investigation.  The  selected  cells  represent  a  set  ot  lines 
detected  from  the  image. 

For  each  detected  line,  a  collection  ot  points  that 
lie  on  the  line  is  gathered  from  the  edge  points,  do 
represent  the  location  of  a  line  in  the  overall  image, 
the  mean  of  the  collected  points,  i.e.,  the  mean  (x,  y) 
value,  is  used.  For  cases  in  which  several  line  segments 
from  different  image  locations  happen  to  be  collinear, 
the  line  segments  will  be  represented  by  just  single 
line  values  of  h  and  0.  In  such  cases,  additional  lines — 
with  the  same  line  values  h  and  0  but  with  different 
mean  {x,y)  values — are  added  to  the  line  set  to  repre¬ 
sent  the  different  line  segments.  The  line  set  is  exam¬ 
ined  further  to  eliminate  those  lines  which  have  mean 
(x,  y)  values  and  h  and  0  values  that  are  close  to  other 
lines.  The  line  set  is  then  used  as  a  set  of  seeds  for 
growing  connected  line  segments  with  a  relaxation 
process,  as  shown  in  Figure  4. 

During  each  iteration  of  the  relaxation  process, 
each  edge  point  in  the  image  is  classified  in  associa¬ 
tion  with  one  of  the  detected  lines.  The  edge  point  is 
then  added  to  the  point  collection  of  the  associated 
line.  When  a  point  is  not  close  enough  to  any  of  the 
lines  (the  allowable  orthogonal  distance  for  point  clas¬ 
sification  should  be  less  than  half  the  distance  of  the 
closest  distinctive  feature  lines  in  the  image)  or  when 
the  edge  direction  does  not  agree  with  the  line  direc¬ 
tion  (the  allowable  angular  deviation  should  be  less 


216  iHf  UMCOU  I'BORSTORV  JOURNll  VOUlMt  6  NUMBIR  1  1993 


•  HSU 

l\\Ihi(.  tin^  / .ifXt  t  /  t  .iiitf  t  s  fittni  tiiul  H^tuXi'  I  ^uppli  i  Iti/axc  ' 


0  (rad) 


FIGURE  3.  Hough  transform  for  straight  lines:  (a)  line  defined  in  x-y  coordinates  and  (b)  Hough  transform  in  h-H  space.  A 
line  in  two-dimensional  (2-D)  space  can  be  represented  by  the  directed  orthogonal  distance  h  to  the  origin  and  the  angle  H 
that  h  forms  with  thex-axis,  as  shown  in  part  a.  For  any  given  point,  the  set  of  lines  passing  through  that  point  can  be  plotted 
in  h-t)  space.  The  result  will  be  a  sinusoidal  curve,  as  shown  in  part  b.  Note  that  the  h  curves  of  two  points  in  x-y  space  will 
intersect  periodically  in  h-d  space.  The  points  of  intersection  (spaced  .r  radians  apart)  correspond  to  the  line  defined  by  the 
two  points  in  x-y  space. 


than  half  the  smallest  angle  formed  by  intersecting 
lines  in  the  image),  the  point  is  not  classified.  At  the 
end  of  each  iteration,  the  current  classification  is  com¬ 
pared  with  the  results  of  the  previous  iteration.  If  any 
difference  exists,  the  line  set  is  recalculated  from  the 
current  point  collection  and,  with  the  updated  line 
set,  the  point  classification  is  reiterated.  The  relax¬ 
ation  process  stops  when  the  current  classification  has 
not  changed  from  the  previous  iteration.  After  the 
relaxation  process  ceases,  the  lines  are  identified  from 
the  edge  points,  and  the  position  and  orientation  of 
each  line  are  calculated  from  the  point  collection  with 
a  least-squares-error  fit.  This  information  is  then  used 
for  simple  shape  formation:  parallel  lines  are  identi¬ 
fied  for  cylinders,  and  intersecting  lines  are  chosen  for 
cones.  Finally,  the  axis  of  symmetry  is  determined  for 
the  position  and  projected  orientation  of  each  chosen 
shape. 

Curve  Extraction 

The  circula’-  base  of  a  cylinder  or  cone  generally  ap¬ 
pears  elliptical  in  angle-angle  imagery  because  the 
acnecr  angle  beween  me  object  body  axis  and  the 


receiver  LOS  is  usually  nonzero.  From  geometry,  the 
ratio  of  the  minor  and  major  axes  of  the  projected 
ellipse  is  the  cosine  of  the  aspect  angle.  The  aspect 
angle,  together  with  the  projected  body- axis  orienta¬ 
tion,  can  be  used  to  determine  the  body  attitude  with 
respect  to  the  sensor  in  3-D  space.  In  extracting  the 
elliptical  base,  it  is  assumed  that  the  straight  bound¬ 
aries  have  already  been  identified.  Thus  the  ellipse 
center  will  lie  on  the  axis  of  symmetry  of  the  shape. 
For  cylinders,  the  major  axis  of  the  base  ellipse  will  be 
half  the  distance  between  the  two  parallel  edges  that 
define  the  shape.  For  cones,  the  major  axis  will  be  the 
distance  from  one  conical  edge  to  the  ellipse  center. 
Thus  the  ellipse-extraction  scheme  first  determines 
the  base  center  along  the  axis  of  symmetry  of  the 
shape  and  then  fits  an  elliptical  curve  at  the  base 
region  to  estimate  the  minor  axis. 

Figure  5(a)  depicts  a  cone  shape  with  an  elliptical 
base.  The  ellipse  function  with  a  coordinate  rotation 
is  expressed  by 

(X  -  (j  -  Jo)r  , 

-  — -  4-  -  ^  ,  (11 
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FIGURE  4.  Relaxation  process  for  line  identification. 

where  (xq,  y^)  represents  the  ellipse  center  on  the  axis 
of  symmetry  having  orientation  and  a  and  h  are 
the  minor  and  major  axes  of  the  ellipse,  respectively. 
Equation  1  can  be  rewritten  for  the  minor  axis  a: 

a  =  ~  ■yo)r| 

,  _(y  -  Jo)? 

■ 

The  above  equation  indicates  that  the  minor  axis  of 
the  ellipse  is  a  function  of  the  location  of  the  ellipse 
center  (xq,  y^)  along  the  axis  of  symmetry,  as  shown 
in  Figure  5(b).  Thus  the  Hough  transform  space 
[a,  {x^|,  Jo)]  can  be  used  for  ellipse  detection  in  which 


values  of  a  are  plotted  for  all  edge  pttints  in  the  base 
region.  Fhe  maximal',  accumulated  cell  in  Hough 
transform  space  will  then  determine  ti  and  the  loca¬ 
tion  of  (.x;,, /„). 

In  addition,  the  edge  detection  of  a  point  ttn  the 
ellipse  curve  should  be  consistent  with  the  curve  tan¬ 
gent  derived  by  differentiating  Liquation  1 ,  as  follows: 

_  ^/(x-X||)^  X  {y  -  y^,), 

tan  f7  =  -  = - ^ - . 

d{y  -  Vii),  b-  (x  - 

Thus  each  edge  point  is  checked  with  respect  to  loca¬ 
tion  and  edge  orientation  before  being  added  to  the 
collection  of  points  for  an  extracted  curve.  The  mi¬ 
nor-axis  value  is  then  refined  from  the  best  curve  fit  of 
the  collected  points.  The  results  of  line  and  curve 
identification  are  shown  and  discussed  in  the  follow¬ 
ing  section. 

Feature  Extraction  for  Angle-Angle  Images 

For  angle-angle  imagery,  targets  are  projected  from 
3-D  space  to  a  plane  perpendicular  to  the  sensor 
LOS,  as  shown  in  Figure  6.  Thus  spheres  are  pro¬ 
jected  as  circles,  and  the  circular  bases  of  cones  and 
cylinders  become  ellipses.  In  the  images,  the  physical 
radii  of  the  circular  forms  of  targets  are  generally 
preserved  without  transformations:  the  major  and 
minor  axes  of  the  elliptical  projection  are  respectively 
the  circular  radius  itself  and  the  radius  with  a  cos« 
factor,  a  being  the  aspect  angle  (i.e.,  the  angle  be¬ 
tween  the  target  body  axis  and  the  sensor  LOS).  The 
boundary  lines  of  targets  are  projected  into  lines  in 
images,  but  their  dimensions  are  generally  transformed 
with  a  factor  of  sina.  With  the  use  of  such  projection 
relationships  between  target  3-D  space  and  image 
2-D  space,  the  size,  shape,  and  orientation  of  a  target 
can  be  inferred  from  its  image  features. 

Two  examples  of  simulated  angle- angle  images  have 
been  analyzed  to  demonstrate  the  extraction  of  target 
features  from  such  imagery.  First,  a  simulated  image 
of  a  complex  object  is  used  for  size  and  shape  estima¬ 
tion.  Then  a  sequence  of  cone  images  is  employed  for 
motion  extraction. 

Size  and  Shape  Estimation 

The  model  used  in  the  example  has  a  cylindrical  main 
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FIGURE  5.  Elliptical  Hough  transform:  (a)  cone  shape  with  elliptical  base  in  x-y  coordinates  after  a  coordinate  rotation  /(,, 
and  (b)  Hough  transform  in  which  a,  the  minor  axis  of  the  ellipse,  is  plotted  on  the  vertical  axis  and  the  location  of  the  ellipse 
center  (xg,  yg)  along  the  axis  of  symmetry  is  plotted  on  the  horizontal  axis.  For  a  detailed  description  of  the  mathematics 
involved,  see  the  main  text. 


body  with  three  cone.s  mounted  at  the  forward  end. 
Figure  7(a)  .shows  a  simulated  angle-angle  image  of 
the  object  at  a  30°  aspect  angle,  and  Figures  7(b),  (c), 
and  (d)  show  the  results  of  a  scene  analysis  that  was 
performed  with  pixel  sizes  of  2,  8,  and  20  cm,  respec¬ 
tively.  For  each  of  the  scene-analysis  images,  the  edge 
points  have  been  detected  and  lines  identified.  (Note: 

1  he  different  lines  are  coded  in  different  colors.) 
Ideally,  two  major  parallel  lines  would  be  selected  for 
the  cylindrical  body,  and  three  pairs  of  intersecting 
lines  would  be  chosen  for  the  cones.  Then  an  ellipti¬ 
cal  curve  could  be  fitted  for  the  cylindrical  base  to 
determine  the  aspect  angle.  Because  the  cylindrical 
body  is  large  compared  to  all  three  of  the  pixel  sizes 
used,  the  rwo  parallel  lines  representing  the  shape  are 
easily  discernible  in  Figures  7(b),  (c),  and  (d).  But  the 
cones,  because  of  their  smaller  size,  appear  distorted 
in  the  images,  particularly  in  Figure  7(d).  Nonethe¬ 
less,  the  cones  are  recognizable  at  pixel  sizes  of  2  and 
8  cm. 

Figure  8  shows  the  quantitative  results.  For  all  pixel 
sizes,  the  estimated  cylinder  radius  (Figure  8[a])  agrees 
well  with  the  model.  F.stimations  for  the  cone  radius 
(Figure  8[b])  and  projected  cone  angle  (Figure  8[c]) 
are  good  for  pixels  smaller  than  8  cm,  and  the  aspect 


FIGURE  6.  Geometry  of  angle-angle  imagery  of  a  cone. 
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FIGURE?.  Size  and  shape  extraction  for  a  simulated  object  consisting  of  a  cylindrical  main  body 
with  three  cones  mounted  at  the  forward  end:  (a)  angle-angle  input  image  after  edge  detection, 

(b)  line  classification  at  a  pixel  size  of  2  cm,  (c)  line  classification  at  a  pixel  size  of  8  cm,  and 
(d)  line  classification  at  a  pixel  size  of  20  cm. 

angle  obtained  from  the  elliptical-curve  fit  (Figure  ratio  for  the  target  at  the  receiver. 

8[dl)  is  also  in  agreement  with  the  model  for  pixels 

smaller  than  8  cm.  Extraction 

In  general,  curve  fitting  requires  finer  image  qual-  Target  motion  can  be  used  to  aid  the  target-classifica- 

ity  than  tine  fitting.  Lines  can  usually  be  detected  tion  process,  as  has  been  demonstrated  in  range-Dop- 

within  one  pixel  to  the  true  edge.  The  pixel  size,  pier  imagery  from  millimeter-wave  (MMW)  and  other 

which  represents  the  sampling  size  at  the  image  focal  microwave  radars.  With  a  time  series  of  target  atti- 

plane,  is  determined  by  the  type  of  focal  plane  and  tudes  extracted  from  such  imagery,  target  motion  can 

the  angular  resolution,  the  diffraction  or  resolution  be  inferred.  The  observation  generally  is  more  straight- 

limit  of  the  imaging  .system,  and  the  signal-to-noise  forward  in  a  time  series  of  high-re.solurion  images. 
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FIGURE  8.  Feature-extraction  performance  at  various  pixel  sizes  for  the  simulated  model  of 
Figure  7:  (a)  calculated  base  radius  of  the  model,  (b)  calculated  base  radius  of  a  cone,  (c)  calculated 
half-cone  angle  of  a  cone,  and  (d)  calculated  aspect  angle  of  the  model.  (Note:  For  comparison,  the 
input  values  used  in  the  simulation  are  indicated  with  straight  lines.) 


For  purposes  of  demonstration,  a  free-body  spin- 
precession  motion  is  considered.  Figure  9  shows  a 
coordinate  system  with  a  spinning  target  and  the 
body  axis  of  the  target  processing  around  the  z-axis. 
The  parameters  0^,  and  p  represent  the  precession 
half-cone  angle,  the  processing  azimuthal  angle,  and 
the  angle  between  the  LOS  and  the  z-axis,  respec¬ 
tively.  Suppose  that  p  and  1)^  are  constant  during  an 
observation.  Then  the  aspect  angle  a  between  the 
body  axis  and  the  LOS  will  be  a  function  of  (p^.  For  p 
in  the  z-x  plane, 

cos  a  =  sin  p  sin  0^  cos  ip^  +  cos  p  cos  .  (2) 

If  ^os«  can  be  estimated  for  different  values  of  (p^, 
then  Equation  2  can  be  used  to  solve  for  the  coefFi- 
cients  sinp  sinff^  and  cospcosffp.  The  sum  and  dif¬ 
ference  of  these  coefficients  are  cos(p  -  0^)  and 
cos(p  +  ffp),  respectively,  and,  from  these  two  quanti¬ 
ties,  p  and  ffp  can  be  obtained. 

The  target  projected  orientation  Oq  in  angle-angle 
images  is  also  a  function  of 


FIGURE  9.  Geometry  of  a  target  that  is  spinning  (a'  out  the 
target  body  axis)  and  processing  (about  the  z-axi  )  simul¬ 
taneously.  The  parameters  Hp  and  (pp  represent  the  preces¬ 
sion  half-cone  angle  and  the  processing  azimuthal  angle, 
respectively,  and  a  and  p  represent  the  aspect  angle  and 
mean  aspect  angle,  respectively. 
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FIGURE  10.  Feature  extraction  for  angle-angle  Images  of  a  simulated  spinning  triconic  body  undergoing  processional 
motion:  (a)  input  images,  (b)  edge  extraction,  (c)  fitted  data  for  cosa,  and  (d)  fitted  data  for  tan  Wq.  The  rotation  rate  chosen 
was  SO'  /sec,  11^  was  16°,  the  precession  rate  for  0p  was  10°/sec,  and  was  90  ’.  The  images  were  simulated  with  a  2-cm  pixel 
size  and  a  2-sec  frame  time  (i.e.,  a  2-sec  interval  between  frames),  with  19  frames  generated  over  a  complete  processional 
cycle. 


FIGURE  11.  Geometry  of  range-Doppler  radar  imagery  of  a  spinning  cone. 
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At  broadside  \  iewin^  (/)  =  I'.cjiiation  }  reduces  to 

tantt,,  -  -  tant^j,  simp^,.  Similar  to  the  case  ot  aspect- 
anyle  data,  the  rtrientatioir  data  can  also  he  fitted  with 
f  iliiation  ,s  tt)  determine  p  and 

As  .m  example,  a  series  ol  simulated  atigle-angle 
images  has  been  generated  tor  a  spinning  triconic 
hodv  undergoing  precessional  motion  (figure  l()|a]). 

1  he  rotation  rate  chosen  was  8()°/sec,  (1^^  was  16°,  the 
precession  rate  tor  0^^  was  IO°/sec,  and  p  was  60°. 
(Note:  The  1,(')S  was  assumed  to  h<-  outside  the  pre¬ 
cession  cone.)  1  he  images  were  simulated  widi  a 
2-cm  pixel  size  and  a  2-sec  trame  time  (i.e.,  a  2-sec 
interval  between  trames),  with  19  trames  generated 
over  a  complete  prece.ssional  cycle,  t  he  triconic  bound¬ 
aries  were  identifted  by  computer  tor  each  image 
trame  (Figure  l()[b|)  so  th.it  the  aspect  angle  «  and 
orientation  could  be  determined.  Fitting  the.se  angles 
to  Fiquations  2  and  5  then  allows  p  and  to  be 
obtained.  For  the  a.spect-angle  history,  Figure  10(c) 
shows  that  the  best  fit  occurs  tor  p  =  91  i°  and 
(1^  =  1 5.4°.  For  the  image-plane  bod}'-a:ds  orientation 
(given  p  =  90°).  Figure  10(d)  shows  that  the  best  fit 
occurs  for  9^,  =  16°. 

Feature  Extraction  for  Range-Doppler  Images 

For  range-Doppler  radars,  targets  with  angular  dy¬ 
namics  are  imaged  along  the  radar  LOS  in  the  D' o 
dimensions  ot  range  and  Doppler.  The  projection  of 
L,  the  length  ot  a  target  measured  along  its  body  axis, 
onto  the  l.OS  is  Lcosa  in  the  rimge  dimension, 
with  ((  being  the  aspect  angle,  as  shown  in  Figure  1 1. 
In  the  ca.se  ot  a  rotating  target,  the  rotational  angular 
velocity  ^projected  onto  the  LOS  is  measured  in  the 
Doppler  dimension  as  Ksino  [7].  With  the.se  rela¬ 
tionships,  the  dimensions  of  a  target  can  be  derived 
from  its  image  range  and  Doppler  extents. 

Application  of  this  technique  is  demonstrated  in 
this  section  with  an  example  of  simulated  range-Dop¬ 
pler  images  of  a  diffuse  cone.  In  the  example,  image 
proce.ssing  and  scene  analysis  are  applied,  as  has  been 
demonstrated  earlier  for  angle-angle  imagery.  A  target 
line  model  is  then  used  to  guide  the  formation  of  a 
target  line  representation  in  the  imagery.  Farget  di- 
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meiisions  and  onem.uions  .ire  then  detetmineil  from 
the  range  and  floppier  extents  ol  a  seijueiiee  ot  the 
imatres. 

htiiUiiiion  of  largi'i 

For  the  cone-sliaped  target  ot  Figure  11,  the  range 
and  Drrppler  extents  aie  related  to  the  target  aspeet 
angle  <<  and  the  rotation  rate  /,  b\ 

R  (h  =  /.  cos  o  ( ^) 


K  ■  z//'  =  4.7/2/,  sin  (t .  (5) 

where  //and  I  Aire  the  range  and  Doppler  extents 
(number  of  image  pixels),  z/rand  z/nare  the  range  and 
Doppler  cell  sizes,  and  /,  and  //are  the  target  phvsical 
length  and  base  diameter,  respectiveb'.  Combining 
the  two  equations  to  eliminate  functions  of  the  aspect 
angle  gives  the  following  ellipse  equation  for  a  rigid 
target  undergoing  steady  rotation  with  constant  /,  D, 
and/: 


4.7/2/ 


The  major  and  minor  axes  ot  the  above  ellipse  can 
be  estimated  for  a  dataset  of  (R,  V)  pairs  by  the  least- 
squares-error  method.  Fhe  physical  length  L  and  the 
Doppler  velocity  can  thus  be  calculated  b\ 


\i^v'  -  ( /?V’ 


L  =  .'== — ^  dr 

\  rW^  -  k'  /?V“ 


4.7 /2/,  =  ,  -  —  dv. 

\  R'V-  -  R-  R-V- 

For  cases  in  which  the  rotation  Irequency  /  is 
known,  the  base  diameter  D  can  also  be  determined 
from  the  Doppler  velocity.  Fheoretically,  the  above- 
derivation  can  be  performed  tor  any  dataset  with 
more  than  two  pairs  ot  (R,  V)  measurements,  pro¬ 
vided  that  the  data  noi.se  or  measurement  errors  are 
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FIGURE  12.  Simulated  range-Doppler  image  of  a  cone;  (top)  before  process¬ 
ing  and  (bottom)  after  edge  points  have  been  detected  and  lines  fitted.  The 
co'or  bar  indicates  increasing  (from  black  to  white)  radar  cross  sectioi 
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nmch  less  than  the  true  data  diHcrcnce  corresponding 
to  ddferent  aspect  angles.  For  a  large  dataset,  the 
nieasureinents  ot  range  and  Doppler  extents  can  be 
median  filtered  separately  to  remove  noise  (while  still 
retaining  the  temporal  trend  resulting  (rom  changing 
aspect  angles). 

With  both  /.  and  Dj^  derived  from  a  given  data.set, 
the  target  aspect-angle  history  can  be  obtained  from 
either  bciuation  4  or  5,  or  from  the  ratio  o(  the  two 
equations.  It  should  not  matter  which  equation  is 
used  to  obtain  the  aspect  angle  because  both  of  the 
equations  were  used  to  derive  F.quation  6,  which  was 
used  to  estimate  both  /.  and  L)^.  More  precise  results, 
however,  may  be  obtained  by  using  the  equation  cor¬ 
responding  to  better  cell  resolution. 

Motion  Extraction  for  a  Simulated  Cone 

For  purposes  of  demonstration,  range-Doppler  im¬ 
ages  have  been  simulated  for  a  cone  (length  of  1 50  cm 
and  a  base  radius  of  19.7  cm)  having  a  diffuse  target 
surface,  A  total  of  280  images  was  generated  over 
28  sec  with  the  cone  undergoing  spin  (spin  period  of 
4.5  .sec)  and  precession  [8].  Each  image  was  53  x  53 
pixels,  with  a  range  cell  representing  5  cm  and  a 
Doppler  cell  representing  2.93  cm/sec.  Figure  12  (top) 
shows  one  of  the  range-Doppler  image  frames. 

Image  processing,  edge  detection,  and  line  fitting 
were  then  performed  on  the  images.  For  the  simple 
known  target  shapes  in  the  current  application,  the 
line  fitting  was  simplified  with  a  piecewise  linear  fit¬ 
ting  of  the  target  boundaries:  the  edges  at  the  base 
were  collected  and  fitted  linearly,  and  the  cone  was 
defined  by  two  lines  fitting  the  cone  edges  and  the 
base  line.  Range  and  Doppler  extents  of  the  target 
could  then  be  estimated  from  the  target  line  model. 
Figure  12  (bottom)  shows  a  frame  of  the  image  with 
the  edge  points  detected  and  the  lines  fitted.  For  such 
figures,  the  bisection  line  of  the  cone  is  the  target- 
body  centerline,  the  distance  from  the  nose  to  the 
base  center  in  range  is  the  target  range  extent,  and  the 
Doppler  spread  between  the  intersections  of  the  base 
line  and  the  cone’s  left  and  right  sides  is  the  Doppler 
extent. 

Figure  1 3  contains  plots  of  the  extracted  range  and 
Doppler  extents  from  the  image  sequence.  Note  that, 
in  spite  of  the  considerable  amount  of  local  noise 


Time  (sec) 

(a) 


Time  (sec) 

(b) 

FIGURE  13.  The  extracted  (a)  range  and  (b)  Doppler  ex¬ 
tents  from  a  sequence  of  simulated  images  (see  Figure 
12). 

present  in  the  imagery,  the  range  and  Doppler  histo¬ 
ries  show  functions  of  coso  and  sinrt,  respectively,  for 
the  aspect  angle  «,  which  changes  because  of  the 
precession  of  the  target.  Target  dimensions  can  be 
obtained  by  the  elliptical  fitting  of  the  range  and 
Doppler  data,  as  described  earlier  in  the  subsection 
“Estimation  of  Target  Dimensions.”  With  such  tech¬ 
niques,  a  value  of  153.5  cm  was  calculated  for  the 
target  length,  and  17.3  cm  for  the  base  radius.  Note 
that  the  derived  target  dimensions  agree  well  with  the 
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FIGURE  14.  Aspect-angle  history  derived  from  the  Doppler 
extents  (see  Figure  13[b]). 

corresponding  input  model  dimensions  of  150  cm 
and  19.7  cm.  Both  image-resolution-related  noise  and 
target-surface-scattering  noise  could  have  contributed 
to  the  differences  in  the  derived  and  input  target 
dimensions.  The  quantitative  impacts  of  such  noise 
are  currently  undet  evaluation. 

Figure  14  shows  the  aspect-angle  history  that  was 
derived  from  the  Doppler  extents.  A  running  averag¬ 
ing  of  4.5  sec  was  applied  to  smooth  the  aspect-angle 
curve,  which  changes  from  38°  to  47°  with  a  period 
of  about  28  sec.  With  the  assumption  that  the  radar 
I. OS  was  outside  the  target  precession  cone,  the  target 
prece.ssion  half-cone  angle  can  be  calculated  as  4.5°, 
which  is  comparable  to  the  input  angle  of  5°. 

Summary 

The  extraction  of  target  features — shape,  size,  orien¬ 
tation,  and  dynamics — from  both  angle-angle  and 
range-Doppler  images  has  been  demonsttated.  The 
target  shape,  size  and  orientation  can  be  determined 
from  single-frame  angle-angle  images,  and  the  target 
angular  dynamics  can  be  estimated  from  a  time  series 
of  such  images.  For  the  case  of  range-Doppler  imag¬ 
ery,  multiple  frames  of  images  are  required  to  obtain 
multiple  measurements  of  range  and  Doppler  extents 
to  enable  the  derivation  of  the  target  dimensions.  A 
time  series  of  range-Doppler  images  can  also  be  used 
to  estimate  the  target  angular  dynamics. 


For  simple  known  targets,  the  fe.iture-extraction 
process  can  be  simplified  greatly  by  obtaining  the 
target  image  description  with  a  piecewise  lineat  ap¬ 
proximation  of  the  target  boundaries.  For  cases  in 
which  some  size/shape  parameters  such  as  cone  half- 
angle,  length,  and/or  base  diameter  are  available,  the 
target  size  and  aspect  angle  can  be  estimated  from 
data  extracted  for  each  frame  of  range-Doppler  im¬ 
ages.  Convergence  of  such  estimations  over  multiple 
frames  of  images  can  then  be  used  to  evaluate  the 
algorithm  performance  and  to  verify  measurements 
for  the  known  target. 

The  algorithm  is  undergoing  extensive  testing  with 
simulated  and  field  data,  and  both  the  precision  and 
efficiency  of  the  algorithm  are  expected  to  improve 
over  time.  Preparation  for  real-time  implementation 
of  the  range-Doppler  algorithm  is  currently  in  progress. 
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